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1 INTRODUCTION 

We focus on Bessel integrals 

f oo 

du u n+1 K (u) K 



with n > and k > 1 integers and on their relation with zeta numbers. We 
present results from [TJ [2l [3l H] and review [5] for a self contained presenta- 
tion. Continous fractions for C(3), C(2) and ip\{l/2>) —ipi(2/3) are presented, 
some of them related to ([1]), others sustained by numerical PSLQ evidence 
[6]. All these continuous fractions are intimately connected to Apery's con- 
tinous fractions for ("(3) and £(2) irrationality demonstrations [7]. Attempts 
to integrate Bessel integrals starting from a multi-integral representation on 
a finite domain are given. In the process Bessel integrals are shown to be 
periods^. Finally a possible way to address the irrationality of ("(5) is pro- 
posed. 



2 REVIEW OF [5] 

2.1 Quantum Mechanics and Bessel Integrals 

The random magnetic impurity model [TJ describes a quantum particle in 
a plane coupled to a random distribution of Aharonov-Bohm fluxes perpen- 
dicular to the plane. It was introduced having in mind the Integer Quantum 



1 Periods are defined in [8j as "values of absolutely convergent integrals of rational 
functions with rational coefficients over domains in R n given by polynomial inequalities 
with rational coefficients". 



J 



Hall effect. A perturbative expansion of the partition function of the model 
in the coupling constant a (the flux expressed in unit of the quantum of flux) 
for 2 impurities, i.e. at second order in the impurity density p, lead us to 
consider Feynmann diagrams at order p 2 a 4 [1] and p 2 a 6 [3]. 

On the one hand, these Feymann diagrams, which reduce after momenta 
integrations to multiple integrals on intermediate temperatures, were shown 
to rewrite in terms of simple [2] and double nested integrals [3] on products 
of modified Bessel functions 



Ip 2 a 4 = / uKq(u) (uKi(u)) du 
Jo 

/*CO f'U 

I p 2 a& =8 / duuKoiufiuK^u)) 2 dxxKi 
Jo Jo 

roo poo 

-4/ duuK^iuKxiu^iuK^I^-uh^K^u)) / dx x K (x) 2 K^xf 

Jo Ju 

POD 

+ / uKoiu^iuK^fdu (2) 
Jo 

On the other hand, one could show |2| by direct integration 

roo q3 i 

J uK a (u) A du = ^— C(3) (3) 
and by integration by part that 

uKoiufiuK^fdu (4) 



is a linear combination with rational cefficienta^] of 1 and J °° uKo^u^du, i.e. 
of 1 and (2 3 — 1)^(3). Similarly one obtained OH] by direct integration!] 

r-oo o 2 — 1 

/ duuI (u)K (uf = —— C(2) (5) 
Jo ° 

and by integration by part that 

duuK {u){uKi{u)) 2 I Q {u) (6) 



and 

poo 

duuKv(u) 2 uKi((u)uIi{u) (7) 



2 Linear combination with rational coefficients means here that there exist three posi- 
tive or negative integers a, b and c such that 



uKo(u) (uKi(u)) du + b uKq(u) du + c = 
Jo 

3 See Appendix A for the derivations of © and ©. 
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are linear combinations of 1 and J °° duuIo(u)Ko(u) 3 , i.e. of 1 and (2 2 — 
1)£(2). Likewise, by integration by part [5] 

oo q roo i roo 

uK {uf(uK l {u)) 2 du = — / uK Q {ufdu--i u 3 K {uf du 
o 15 Jo 5 Jo 

(8) 

It was then natural to argue [5] that ([2]) might also rewrite as a linear com- 
bination with rational coefficients of simple integrals on product of Bessel 
functions of weight 6 -defined as the total power of Bessel functions- and of 
zeta numbers of weigth 6 — 1 = 5, like £(5), or below. The "counting rule" 
inferred from ([31 HI El El El (HI) is that an integration du u n ( ) with n 
odd diminishes the Bessel weight by one, so f °° u Kq(u) 4 (uKi(u)) 2 is like 
C(5), that a /-Bessel function has no weight, so J Q °° duuIo(u)Ko(u) 3 is like 
C(2) and duu 3 Koiu^K^u) 2 dxxK^h^Koix) 2 like C(3)C(2) he. 
like C(5), that in turn is like J^°uKq(u) 6 or f °° u 3 Ko(u) e -why only odd 
powers u and u 3 appear here and no higher power will become clear later. 
Indeed a numerical PSLQ [6] search gave 

-y roo ^ f°° 2^ 1 

W =PSLQ ^ J u K (uf du +^j o u 3 K {uf du - -j^- C(5) (9) 

(the 2 5 — 1 factor multiplying £(5) has to be viewed in parallel with 2 3 — 1 
multiplying £(3) in ([3]) and 2 2 — 1 multiplying £(2) in ([6])). Note that from 
now on an identity obtained from a numerical PSLQ search will be labelled 
=PSLQ as in ©. 

The fact that the double nested integrals in ([2]) can be reexpressed as linear 
combination with rational coefficients of simple Bessel integrals with the 
right weight -here weight 6- and of £(5), fits well in the "Bessel integral — > 
zeta number" mapping. If a product / of Bessel functions is, like in ([3l HI El 
El El) j mapped by simple integration on a zeta number denoted by £(/) 

oo 



/-> / f(u)du = ((f) 
Jo 

then for a pair of such products /, g the mapping by double "nested" inte- 
gration 

poo ru 

f,9^ f(u)du g(x)dx = ((f,g) 
J J 

on a polyzeta number denoted by Cif id) makes sense, since, because of 

oo ru poo roo roo ru 

f(u)du / g(x)dx = / f(u)du / g{x)dx — / g(u)du / f(x)dx, 
o Jo Jo Jo Jo Jo 

one has 

C(/,s) = C(/)C(s)-C(5,/) (io) 
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in analogy with 



C(p,q) = C(p)C(q)-C(p + q)-C(q,p) 



(11) 



for the standard polyzeta £(p, q) = Xm>m is"^ C(p> <?) would be defined 



as En 



1 1 

>m nP mi 



+ IC(P + q) then (jnj would take the form ([TO]) . 



2.2 Recurrence 

It appears that 1 and J °° u K^u)^ du one the one hand, and 1, J °° u Ko(u) e du 
and Jq 00 v? Kq(u)® du on the other hand, play the role of building blocks in 
the linear combinations reexpressing, for example, Q and ([8]). This might 
indicate their special role as basis for more general families of Bessel inte- 
grald3- 



Indeed, consider such a family of Bessel integrals with a given weight re 

1 



I 



(«) 



n n+i K ('u)' t ~ J ^i(M) J du 



0,1,..., K 



(12) 



where n > re — 1 for 1^ to be finite. Integration by parts gives the mapping 

(13) 



j{ K ) , r(") 
n i K — ^ -'n+l,K 



7-00 



n + 1 



n - jf + 2 
with the (re + 1) x (re + 1) matrix 



(*) 



/ n K ( n+1 ) 

' u n+2 



1 






V o 





2(n+l) 
n 











K - 1 



3(n+l) 
n-1 











(«-2)(n+l) 
n 



4(n+l) 
n-2 










(re-3)(n+l) 
n-1 
















«(n+l) 
n— re+2 



\ 





n+1 



whose determinant when re is odd is 
even is vanishing since in this case 

re/2 



V(-l)<(n-2/ + 2) TU 



(3) 2 (5) 2 ...(K) 2 (n+l) 2 ''+ 1 
(ra+2)(ra+l)(ra)...(ra-re+2) > 



«/2\ r (re) 



1=0 



I 



L n.2l 



n— re+3 

o / 

(14) 

and when re is 



(15) 



4 The same is true of 1 and du uIq{u)Ko(u) 3 for and more generally for the 
family of Bessel integrals obtained from (|12|) by replacing either a Kq by a Jo or a K\ by 
a 7i, see [3]. 
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Applying (|13p twice gives the mapping n — > n + 2 

r W (n + l)(n + 2)(j-l)j (K) 

n '^' (n- j + 2)(n- j + 4) n + 2 J" 2 



(w + l)(n + 2) 
n - j + 2 



n — j + 2 n — j + 4 



i n+2,j 



( w + l)(n + 2)(«-j-l)(K-j) W 
+ (n-j + 2) 2 ™+ 2 ^ +2 1 bj 

(re) 

which conserves the parity of n — j. It follows that the I^l's are divided into 
two sub- families, depending on the parity of n — j. 

Let us focus on the sub-family n — j ever0: it is enough to assume that n is 
even, and consider 

iQ, J$, Jf], n>K (^even) (17) 

or 

4?. 45' J S-l> n>K-l ( K odd) (18) 
and note that 

4Si,i, 4+1,3: 4+l.ic-l («even) (19) 



or 



4Sl,3. 4+l,« (^°dd) (20) 



are respectively related to (|17p and ()18p by inverting (|13]) . keeping in mind 
(fTBD when k is eveiM 



By inverting f)16|) (again keeping in mind (|15p ). all integrals in (|17p (and 
thus in (|19|) ) are linear combinations with rational coefficients of the initial 
conditions 

I ( K % l£>} ( K even) (22) 

Now, on the one hand from (|13p one has 2lj^}_ 1 K _ l = n(lj£l + (k — 1)4k— 2) 

and, on the other hand, „_i = 1/k\. It follows that 1^2 can be replaced 
by 1. Using (|15[) . one can drop one more element: therefore, all integrals in 
(fl~7]) (and in (|19p ) are linear combinations with rational coefficients of the 
k/2 numbers 



An analysis for the sub-family n — j odd can be done along the same lines. 
6 The integrals (JT7J) and (|19p (respectively (|18|) and l|20|)l are tantamount to the set 

> 

u^Koiuf-^uK^ujy du neven>0 (21) 
with n even (respectively k odd). 
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{1, 45 > 41 4:i-4l («even) (23) 

Likewise, all integrals in (j 1 8|) (and thus in (J20j) ) are linear combinations with 
rational coefficients of 

Ul-V ^1,2, -, (^odd) (24) 

that is to say, since I^-i K -i = V^'i °f the (k + l)/2 numbers 

{1, l£! li0) 1^,2, (^°dd) (25) 

Finally, by applying ()13p appropriately for < n < k, ()23p can be mapped 
on 

{1, ig, ig, } (Keven) (26) 



and fl25J) on 

{1, ig, /g, I^ 3i0 } (*odd) (27) 

It follows that for n even the /t/2 numbers in ([26|l . namely 

{1, / ^o(w) K dn, / A„(«) fi , / n K - 3 Ko(«) K } (28) 

io jo Jo 

constitute a basis for the integrals (fT7|) and (|19p . Likewise for k odd, the 
(k + l)/2 numbers in (|27p . namely 

/•OO /'OO /'OO 

{1 , / uK Q {uYdu , / U 3 ^ (n) K ,...,/ u K - 2 K Q (ur} (29) 



constitute a basis for the integrals (|18p and (|20p . 

By basis one means that these numbers should be independent over Q: none 
of them is a linear combination with rational coefficients of the others i.e. 
there are no positive or negative integers a, a n such that 

AC— 4 t*OQ 

a x 1 + a n x / u n+1 ET (u) K = (k even) (30) 

n even=0 ^ 

or such that 

/c— 3 

a x 1 + a n x / n n+1 ET (u) K = (« odd) (31) 

n even=0 

This implies that the numbers in (|28p with k even are irrational and also 
irrational the one relatively to the other. The same should be true of the 
numbers in (1291) with n odd. 
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2.3 Asymptotic eigenvalues 



In the asymptotics limit n — > oo, the mapping (|16p becomes 
45 = (j-l)j/M 2]i _ 2 + (2j^-j^^ (32) 

The [(k— k mod 2)/2 + l] eigenvalues of the resulting matrix are k 2 , (k — 2) 2 , 

(k — 4) 2 , . . . where the last eigenvalue is 1 for odd k and for even k. In the 

latter case one has to reduce the dimension of the matrix by one unit, using 

(fT5|) . Thereupon, the largest eigenvalue of the inverse matrix is 1 or 1/4 for 

k odd or even, respectively, whereas the smallest one is l/n 2 . Remarkably, 

with the initial conditions (|22p . (f24"|) . it is the smallest eigenvalue 1/k 2 that 

determines the asymptotic behavior of the inverse matrbfj This can be 

understood by noting that the eigenvector corresponding to the eigenvalue k 2 

is {1, 1, ... , 1}, since (j— l)xl+(2j(n— j)+K)xl+(n— j— l)(/e— i)xl = k 2 for 
(k) 

all j's. Now, the /„4-oo j' s happen to not depend on j, because their integrand 
u n + 1 Kq{u) k ~ j 'Ki(u) j peaks when n — > oo at large values of u, where both 
Kq(u) and K\{u) are approximated by K v (u) — ->■ \JT^e~ u - Therefore, in the 

asymptotic limit, the vector {/^^ •} is indeed proportional to {1, !,...,!}. 



2.4 Question 

The irrationality claims on (|28p and (|29p have been checked numerically. 
They are also supported by the eigenvalues discussion above: (f22|) and 
(|24|) uniquely lead to an asymptotic behavior governed by the correspond- 
ing smallest eigenvalue and are decomposed on the basis (|28p and (|29p . 
It might mean that the building blocks of (f28|) and (|29p . Bessel integrals 
J °° u n+l Kq{u) k du with n even, play some special role in number theory. 
Clearly the cases k = 1 and k = 2 are trivial since u n+l Kq{u) 1 ~^ {uK\{u)y du 
with j = 0, 1 and J* °° tt n+1 ^(w) 2- - 7 (^-^l (uW du with j = 0, 1, 2 are all ra- 
tional and accordingly both the basis (|27p for k = 1 and (f26|) for k = 2 have 
for sole element 1. So the question asked when k > 3: is it possible to assess 
the irrationality of the basis (|28p and (|29D ? 

One will show in the next section that in the first non trivial cases K = 4 and 
K = 3 where the basis (|26[) and (f27|) have dimension 2 the algebra above leads 
to continous fractions with seemingly insufficient fast numerical convergence 
to hint at the irrationality of /q 4 q or /q 3 ^ . 

As far as the case k = 4 is concerned, this question is formal since Jq 4 q 
can be integrated to a rational number times C(3), which is known to be 

7 Starting from any other initial condition would lead to an asymptotics governed by 
the highest eigenvalue. 
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irrational [7]. When k = 3 on the other hand, Jq q can be also be integrated 
to a rational number times ipi(l/3) — , 0i(2/3), a £(2)-like number, whose 
irrationality is so far not known. Some intimate relation with Apery's proofs 
of the irrationality of £(3) and £(2) will show up in the process. 



3 CONTINUOUS FRACTIONS 



3.1 Weight k = A — y £(3) 



The n — > n + 2 mapping ([16 



/ j(4) 
' 1 n,0 
j(4) 
J n,2 

V I (4) 



2(1 + n) 



(i) 



' 2+n 





6 

2+n 
6(l+n) 
5 — 

n^ 
6(2+n) 





(2+n) 
^ — 

2(2+n) 



leads to, using (n + 2)7^ 



\ u (-2+n)n (-2+n)n 

2nig + (n - 2)7$ = in (JS 



/ j(4) 
/ 1 n+2,0 
j(4) 
-'n+2,2 
I j(4) 
\ I n+2A 




+ 

with asymptotics when ra — )• oo 
6 10 



32+10n 
(2+n)* 
6(4+n) 



6n 
(2+n)* 
8+10n 



2 

A+-A- 



2 

A++A- 



(-2+n)n (-2+n)n 




(33) 











and eigenvalues {A_ = 4, A+ = 16} and eigenvectors {{1, — 1}, {1, 1}}. As 
already stated in the Introduction the asymptotic scaling is governed by the 
asymptotic eigenvector {1,1} corresponding to the eigenvalue A+ = 16 (and 
thus to jq for the inverse iteration). Since n is even set n = 2k and define 



u(k) 



I. 



(4) 
2fc,4 



f (4) 
l 2k,0 



(34) 



to obtain 
u{k) 

with 



+ 



a(k)d(k) - b{k)c(k) 



a(k) + c(k)u(k + 1) _ c(k) 

b(k) + d(k)u(k + 1) ~ dQcj ' d(k)(b(k) + d(k)u{k + 1)) 

a(k) = 3(2 + k)(l + k) 2 
b(k) = (8 + 5fe)(-l + k)k 
c(jfe) = (2 + 5/c)(l + A;) 2 
d(fc) = 3ife(-l + fc)Jfe 



(35) 
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Call 

d[k)u{k) - c(k) = z(k) 
and get the iteration \z{k) — > z(k — 1)] with the continuous fraction 

(a(k)d(k) - b(k)c(k))d(k + 1) 1 



z(k) 



d(k) ^ + l)(S} + §S}) + *(fc + l) 

This is finally 

Z{k ~ 1) = 2 + 9k + 15k 2 + Wk^ + z(k) (36) 
with k — ^ oo asymptotics 

-16/c 6 , z(k) 

v ; 10fc 3 + z(fc) fc->oo k 3 



here lim^oo z(k)/k i = lim^-^oo (d(k)u(k) — c(/c))/A; 3 = —2 since one iter- 
ates from it (oo) = 1. Calling z(k) = d(k)(^P^ + ^^) + z(k) = d(k)(u(k) + 
fe(fc-i) - 



d(fc-i) ))' ^ e mverse iteration is 

-16/fc 6 



5(fc) 



2 - 9/c + 15/c 2 - 10/c 3 + - 1) 
with asymptotics 

^(fc) - 77~rq TTTT =>• hm = {2, 8} 

v ; -10/c 3 + z(fc) fc^oo A; 3 1 ' J 

here lim^oo z(k)/k 3 = lirm : _ > . 00 (d(k)u(k) + b(k))/k ?J = 8 since one iterates 
from the Bessel integrals initial condition z(2). 

Continuous fractions appear in Apery's proof of £(3) (and C(2)) irrationality: 

C(3) = , 6 + °" r (37) 

5 + 4rf^ 

P(2)H d 



P(3) + ... 

with the Apery polynomial 

P(k) = 5 + 27k + 51k 2 + 34k 3 = (2k + l)(17k 2 + 17fc + 5) 
This corresponds to the iteration 



— A; 6 
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and to the rational approximation 



C (3) = lim ^MhMlMl (39 ) 

fc^oo 2/(KJ{j / (o)=1,j/(1)=5} 



with 



y(k + 1) - P{k)y{k) + A: 6 y(A; - 1) = (40) 

The initial conditions are i.e. y(0) = 0, y{\) = 6—7- y(A;)| y ( )=o,y(i)=6} an d 
2/(0) = 1, y(l) = 5 -> y(fc){ 2/ ( )=i, 2/ (i)=5}. C(3) is proven to be irrational 
thanks to the sufficiently fast convergence with respect to the increasing size 
of the denominators of the rational approximation (|39p . 

In the weight k = 4 case the iteration (f36|) starts at k = 2 (n = 4). From 



z{2) = -16 x 3 6 = Hm y(k) {y{2)=1M3)=0} 

2 + 9k + i 5fc 2 + 10 p| fc=3 + 2+9fc+15 -^ u _ 4+ ^ i/(*) {y(3)ss0iy( s)=i} 

(41) 

with 

y(fc + 1) - (2 + 9fc + 15/c 2 + 10fc 3 )j/(fc) + 16/c 6 y(A; - 1) = (42) 

Using ()13p the k = 2 (n = 4) initial conditions are reexpressed on the basis 
([28]) for k = 4, i.e. on {1, 1$}, as 

7 (4) = __L xl + JL x7 W and 7 (4) = J3_ J_ (4) 

4 '° 512 384 °>° 4 ' 4 1536 128 0,0 

so 

z(2) = -96 37 " 36 t° = lim y( ^fa(2)=i,y(3)=o} (43) 
27 - 28^ o yv fc ){j/(2)=o,j/(3)=i} 

It appears numerically that the convergence of the rational approximation 
(|4ip does not seem sufficiently fast to prove the irrationality of z(2) that is 



to say that of /q 4 q . One notes in (|36p 

• in the denominator 2 + 9A; + 15k 2 + 10A; 3 = P(k) - 3(2fc + l) 3 where 
P(k) is the Apery polynomial. 

• in the numerator k & as in the Apery case. 
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These are not coincidences: Jq 4 q = J °° uKq(u) a du has already been inte- 



>0 

grated in ([3]) to be proportionnal to £(3) 



\poo 1 _ v-^oo (— l) p 

45 = £^i <( S) = ^"-^-^^ (44 ) 
so that (|43p is in fact a rational approximation to £(3) 



z(2) = DG 74 " 63 ^ = hm y(fe) faW=i.v(3)=Q} 
54 - 49C(3) fc^oo y(fe) {w( 2 )= o )V(3)=1} 

One can push further the analogy by rewriting both Apery and weight k = 4 
rational approximations to £(3) starting from z(0) i.e. with initial conditions 
{y(0) = 1, y(l) = 0} and {y(0) = 0, y(l) = 1}. 

Apery: from 

= 6 + OQ = Hm f( fc ){y(0)=0, V (l)=6} 



5 + 1 



16 \ k ^°° 2/(fc){y(0)=l,y(l)=5} 



one deduces 

- (0) 6 5 ~ 16 lim y( k ){y(Q)=iMD=o} 

C(3) P(l) + n/nx 7 26 _ 36 fc^oo y(fe) {s ( )=fl lV (i )= i} 

F \ 2 )+ P(3)+... 

with y(fc) given in (|40p . 

Weight k = 4: one starts from z(2) in (j4"3|) and with (|36p iterates to 

*(0) = -4 - 2 (45) 

^0,0 

where a simplification has occured since one would have expected in general 

*(0) 



c + dl$ 



Finally one obtains 



* (0 )=*- 2 = -^-2 = lim yi ^y(o)=i,vm=o } 

2/ ^ ' k ^°° y^){y(0)=0,y(l)=l} 

with y(k) given in (j42]). 
So finally one has 
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Apery: 

— k® —k G 
Z ^ k ~ l > = P(k) + z{k) = (2k + 1)(17A; 2 + 17k + 5) + z{k) 

with asymptotics 



~ oa,i ^ 7Ts => lim 4^ = {-(! + ^ 

v ; 34/c 3 + z(/c) fe->oo k 3 1 v y (l + \/2) 4 



y(fc+l)-P(%(fc)+£; 6 y(fc-l) = y 2 -34y+l = -> {(1+V2) 4 , ^^I > 

,(0) = -^--5 = Um yWfa(0)=i,y(i)=0} 
7C(3) fc^oo y(fc){ 2/ (o)=o, 2 /(i)=i} 

Weight^] « = 4: 



Z( ' ~ - 3(2fc + l) 3 + ~ (2fc + l)(5fc 2 + 5A; + 2) + i(jfe) 

(47) 

with asymptotics 

-16/c 6 , z(k) 

v ; 10fc 3 + z(k) fc^oo A; 3 



y(fc+l)-(P(fc) - 3(2fe + l) 3 )y(k)+16k 6 y(k-l) = -►*_>«, y 2 -10y+16 = -> {8, 2} 



8 One could use that 2 + 9fc + 15k 2 + 10k 3 = P(k) - 3(2fc + l) 3 = (1 + 2k)(5k 2 + 5k + 2) 
is a multiple of 2 : then 

-4k 6 —4k G 

<■ J P( fc )-3( 2fc + l)3 +z(fc) (2fc + l) 5* 3 +5*+2 +g(fc ) 



2 



y(/c + 1) - ^ 3 2 ^ fc + 1 ^ y(fc) + 4fc 6 2 /(fc - 1) = -> y 2 - 5y + 4 = -> {4, 1} 

2(0) = - 1 = lim ,W W» =1 '»'" = °> 

7C(3) fc-^oo y(fc){»(o)=o,»(i)=i} 
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z(0)= 12 _ 2= lim y(k) {y{0)=1M1)=0} 

7C(3) fc^oo y(fe){ 2/ ( )=o, 2 /(i)=i} 

A PSLQ search also gave 

—A; 6 —A; 6 
^* ~ = P(k)-2(2k+i)3 = (2k + l){Zk 2 + Zk + 1) + z(k) 

with asymptotics 



^(fc+l)- P(fc) 2 ( 2fc + 1 ) y {k)+k*y(k-l) = y 2 -Qy+l = -> {(l+^2) 2 , ^} 

3 (1 + V2) 2 

*((\\ - 8 1 - Hm ^ fc W)=i,y(i)=o} 

z[v) -pslq 7777^- - 1 - lim — 

(C,{6) fc^OO y{K){ y (p)=o t y(i) = i} 

(one has used that P(k) - 2(2k + l) 3 = (2k + l)(9k 2 + 9 A; + 3) is a multiple 
of 3). 

Only the Apery rational approximation has a sufficiently fast convergence to 
check numerically the irrationality of z(0) = — 5, i.e. of C(3). 

3.2 Weight k = 3 ^i(l/3) - ^i(2/3) 



The n — )• n + 2 mapping (|16p is 

/ 7 (3) \ / _A_ _6_ \ / j(3) \ 

J n,0 _C-| , „\ 2+n 2+n 1 n+2,0 
A3) I — (1 + nj I 2 (6+7n) (3) 

\ ^71,2 / \ n / \ ^n+2,2 / 

with asymptotics matrix n — ?■ oo 

3 6\_ / *±±^= 3(a + -a_) \ _ / 1 1 \f A+ \f 1 1 
2 7 J " ^ ^ ^ J"^ _i ) { o A J V 1 -i 

with eigenvalues {A_ = 1, A+ = 9} and eigenvectors {{1, —1/3}, {1, 1}}. As 
already stated in the Introduction the asymptotic scaling is governed by the 
asymptotic eigenvector {1,1} corresponding to the eigenvalue A+ = 9 (and 
thus to | for the inverse iteration). 
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Define 



call 



u{k) 



r(3) 
1 2k,2 

1 2k,0 



z(k) = d(k)u(k) - c(k) 



and obtain, using 



a(k) = 2k(l + k) 
b(k) = 3k 2 
c(k) = (l + k)(3 + 7k) 
d(k) = 6k 2 

the iteration [z(k) — > z(k — 1)] with the continuous fraction 

-9A: 4 

Z{k ~ 1)= 10k 2 + 10k + 3 + z(k) (48) 

with asymptotics 

/ \ -9/c 4 , z(k) 

z(k) =i — ro 77T im 4/ = {-1, -9} 

v ; 10k 2 + z(k) fc-^oo /c 2 i ' J 

here — /c 2 since one iterates from u(oo) = 1. Calling z(A;) = 2 + 9k + 15/c 2 + 
10A; 3 + z(k), the inverse iteration [z(k — 1) — )■ is 

z(k) 



-3 + Wk-Wk 2 + z(k-l) 
with asymptotics 

z(fc) " -io fc 2 + ^) ^^l^ = {1 ' 9} 

here 9/c 2 since one iterates from the Bessel integrals initial condition z{l). 

Using (|13p the k = 1 (n = 2) initial conditions are reexpressed on the basis 
(f29j) for k = 3, i.e. on {1, I^o}, as 

1 9 i f (3) 1 

I (3) - -- x 1 + - x / {3) / (3) - - x 1 => u(l) = 3i - 1 (ao) 

1 2,0 zi "^-'0,0 



*(1) 



-23 + 40/, 



(3) 

0.0 



-1 + 2/, 



(3) 



z(0) 



-3 + 



0,0 



2/ 



(3) 
0,0 
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where again a simplification has occured as in (j45J). 
So the rational approximation 

<0 ) = -3 + * lim ggjgfcMzS (50) 

2I o,o y( fc )b(o)=o,^i)=i} 

with 

y(fc + l)-(10/c 2 + 10/c + 3)y(/c) + 9A; 4 y(A;-l) = -> y 2 -10y + 9 = -> {9, 1} 

(51) 

The polynomial 10/c 2 + 10/c+3 = llfc 2 + ll/c+3— fc(fc+l) in the denominator of 
(|48|) as well as /c 4 in the numerator (with the opposite sign) are again "Apery 
like". Indeed Apery's proof of £(2) irrationality reads 

z(k-l) 



z(Q) 



11/c 2 + 11/e + 3 + z(/c) 
5 



C(2) 



(3) 

This is not a coincidence: for weight k = 3 it is possible to integrate Iq q 
Jq 00 uKq(u) s du to be proportionnal to a £(2)-like number^ 

^(1/3) -Vi (2/3) 



uKo(u) 6 du 
keeping in mind that 

c(2) = ^(1/3) + ^(2/3) 
where tp n {z) is the PolyGamma function 

oo 1 

^(.) = (-ir +1 n!^ (52) 



Again the numerical convergence of (|50p does not seem sufficiently fast to 

{ 3) 

conclude to the irrationality of z(0), that is of Tqq, that is of ■0i(l/3) — 
ipi(2/3), whose irrationality is unknown. 



There are two other cases (see van der Poorten's review [7]) 



<k - 1) = n — 2(0) = 7t?jy — 2 

v ; (7k 2 + 7k + 2) + z(k) ( ' c(2) 

and 

z(k p.- +fc 4 (4fc + l)(4fc-l) ^./rn_ 5 _ l 

Z{k ij - (2k + l)(3k 2 +3k + l) + z(k) ^ ZW ~ 1 

We found another case of the last type by a PSLQ search, see Section (|3.3[) below. 
10 See Appendix A. 
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3.3 Continuous fractions summary 



• C(2): Apery, two other cases in Van der Poorten, PSLQ nu- 
merical 

z(k -D- +fe 4 (4fc + l)(4/c-l) _ _j_ _ 

Z[k ij " (2k + l)(3k 2 + 3k + 1) + z(k) ^ ZW ~ 2 C( 2 ) 

n t\ +3k 4 (3k + l)(3fc - 1) , nS 7 

Z{k ~ l)= (2k + l)(13k 2 + l3k + 4) + z(k) Z{0) =PSLQ " 4 

. Vi (1/3)-^! (2/3) :k = 3 

Z{k ~ 1)= (Wk 2 + Wk + 3) + z(k) ^ Z{0) = ^i(i/3) 1 4(2/3) " 3 

• C(3): Apery , PSLQ numerical, k = 4 

" 1} = (2fc + i)(i7^ + m + 5) + z(fc) ^ z(0) = c&y " 5 
" 1} = (2 fe + i)(3^; / 3l + i) + ,( fc ) z(0) =PSLQ - 1 
z(fc ~ 1) = - ^pfS : =>^) = ^y-i 

(2fc + l)( + l)+*(fc) 

where one notes 

• the usual pk 2 + pk + q in the denominators and q in the z(0)'s 

• the sign changj^l +/c 4 — > —k A in relation with ipi(l/3) + ?/>i(2/3) — > 
^i(l/3)-Vi(2/3). 

Various PSLQ searches did not produce so far any other example of continu- 
ous fractions of this type for £(2), C(3) or ■0i(l/3) — -01 (2/3) (for which only 
the casd^H listed above seems to be known) . 



11 See also (|44[): from this point of view Apery's - u uu U m i^ul^l i 0n u ro <a,o 7| ,^ 
12 This continuous fraction appears as a "sporadic" case in the numerical search [9] where 

it is related to the Dirichlet L-series £3(2) which "is, up to a factor 3 3 ^ 2 /4, the maximum 

volume of a tetrahedron in hyperbolic 3-space". 
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4 TRYING TO INTEGRATE f™ du uK {u) n 

As alluded to in Section (J3J) Feyman diagrams momenta integrations lead to multiple 
integrals on differences of consecutive intermediate temperatures. Bessel integrals 
J duuKo(u) n can quite generally be represented in this particular way: change 
variable u = 2\ft and use the integral representation 

daa^e^-i = 2K v {u){^) v (53) 



so that 



du uK (u) n 



n 



OO />00 l />oo i -i 

/ den / da n 1 - ai --L -o 2 -7f 
2dt / — — . . . / — e a i e a ^ . . . e 



o Jo 



2 J 2 a\a 2 ... a 



/* OO 

=2 



dai Z" 00 da„ 1 1 



o 2 "Vo 2 ai a 2 ...a„ i + JL + . . . + JL 

■1 /*00 /*00 -j^ 



ai a2 

OO 

(ai+a 2 +...+a n ) 



dai ■ • ■ / da n 

2 n 1 ,/ J • ■ • a,n-i + a 2 a 3 • • • a n + . . . + a n ai . . . a n -2 



(54) 



a n integration: introduce the variable /3 to rewrite e ("i+ a 2+-+»J — d/3 e ^S(P~ 



(ai + «2 + • • • + an)), then integrate over 
duuiT (u) n 



X 



2 



n-l 



/./3 rfi—ai rf)— ai— o„_2 

d/Je - ' 3 / dai / da 2 • • • / da„_i 
Jo Jo Jo 

1 



aia 2 . . . a„_i + (f3 - (ai + a 2 + . . . + a n _i))(a 2 o 3 . . . a„_i + 0,3(14 . . . a n -iai + • • • + a>\ ■ ■ ■ ^-2) 

(55) 

Change variables — > a' t = at/ ft (notation a' t — > a^) so that the /3 integration 
becomes trivial and finally 

00 ^ pi pl — ai pl — ai—...—a rl — 2 



duuK (u) n = -^—f I dai / do 2 . . . / da n _i 



2 jo 

1 



aia 2 ■ ■ ■ a„_i + (1 - (a± + a 2 + . . . + a„_i))(a 2 a 3 . . . a„_i + a 3 a 4 . . . o n _iffli + . . . + aia 2 . . . a„_ 2 ) 

(56) 

With the notations 
u„ = ai + a 2 + . . . + a„ 

t>„ = a 2 a 3 . . . a„ + a 3 a 4 . . . a n a\ + . . . + a„aia 2 ■ ■ . a„_ 2 + aia 2 . . . a n -i 
u>„ = aia 2 . . . a„ 
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one has shown that 



duuK (u) n = — - / doi / do 2 ... / da„ — e u " 

^ JO JO JO v n 

1 /"^ rX—Ui pX—u n —2 



= T^ZT / dai / da 2 . . . / da n _i 

2" 'Jo Jo JO Wn-l + (1 - U n -l)«n-l 

(57) 

a„ i integration: use 

«n-l = a n-l + Un-2) «n-l = 0„_1?;„_2 + U>„_ 2 , U>n-1 = 0„_lW„_2 

so that 

1 1 



U>„-1 + (1 - «r»-l)«n-l w«-2(l - u„- 2 ) + a„_iu„_ 2 (l - Un-2) - a^_ 1 w„_2 

1 1,11 



-<u n _ 2 (a„_i - a+_ 1 )(a„_i - a n _ x ) -u„- 2 (a+_i - a„_i) a n _i - a+_ x a„-i 

(58) 



where 



-v n - 2 (l - u„_ 2 ) ± \Jvl_ 2 (l - Un-2) 2 + 4w n - 2 v n -2(l - Wn-2) 



One finds 



1 



1 p 1 — u\ rX—u 



On _ 

z JO JO 



duuK (u) n = — / dai / da 2 . . . / da„_ 2 



1 1 I a n-l 1 ~ U «- 2 ~ a n-l 1 

log |-+ — -^| (59) 



^Jvl_ 2 (l - U„_ 2 ) 2 + 4w„_ 2 Vn-2(l - W„- 2 ) a »-! 1 Un - 2 

where the integrand rewrites as 



1 i a »-i| 2 1 + X 

log I — r — I = 7- ; Xlog- 



v 2 n _ 2 (l-u n - 2 ) 2 +4w n - 2 v n -2(l-u n ^) <-i (1 " «n-aK-a 1 



1 



(60) 

with 

*-./= 



(1— U„_ 2 )t)n-2 



Change variables a% — > Xi — Xai/u n - 2 so that when < a,; < 1 then 1 > X > 
1 > x; > 0. The Jacobian is 



(u n - 2 ) 



n-2 
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Use homogeneity relations 

n— 2 



i=l 
n-2 



^ aidiV n -2 = (n - 3)u„_2 
i=i 

n-2 

^ a l 9iW„_ 2 = (n - 2)u;„_ 2 (61) 



so that 



and 



Eai = -ki-i 2 )-^ 

r-f 2 1 - u„_2 



/ dwwA'o(u)" =t— f / dxi / dx 2 ... / dx„_ 2 log- 

Jo ^ Jo Jo Jo 



/o Jo Jo 

2 K- 2 )"- 2 



1 + 1 

1 -X 



(i- u „_ 2 K_ 2 x»-3|-ix(i-x2) T3 J- ; | 

=^rr / da; i / ax 2 ... I dx„_ 2 log 



o Jo 



1-Xv n _ 2 X"-3(l-X 2 ) 

(62) 



Use 



Xn-2 / \n— 2 

= K-2j xix 2 ...x„_ 2 

J_ _ 4w»_ 2 ,„„., 

X 2 (1 - un-ajwn-a 1 ' 



so that 



f°, i r 1 , f l - xi A r ' • . l+n i-w n _ 2 

/ ditiiAou =— - / dxi / dx 2 ... / dx n _ 2 log- - — 

Jo 2 Jo Jo Jo 1- X X XiX 2 ■ ■ ■ X„_ 2 



Use 



X =X\ + X 2 + • • • + X„_2 
U„_2«n-2 1,1 1 



W n - 2 X X! X 2 X n -2 

and (|63|) rewritten as 



-X 



(1 - U n _ 2 )w„-2 (x 1 +X 2 + ... + X„_ 2 ) 2 

so that finally 

°° , 1 Z- 1 -* 1 , /• ' ' ' , l+X 1+ X 2 + ...+X rl _ 2 

duuK a (u) =7-—r dxi dx 2 ... / dx n _ 2 log ■ 



^ Jo Jo Jo 

4 



1 - (xi + X 2 + . . . + X„_ 2 ) 



4(xi + . . . + x„_ 2 )xi . . . X„_2 + (1 - (xi + . . . + x„_ 2 ) 2 )(a: 2 2:3 • ■ • X n -2 + • • ■ + XiX 2 . . . x„_ 3 ) 

(64) 
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This integration generalises to du u v K§(u) n with p odd: for example 

duu 3 K (u) n = / dai / da 2 . . . / da„_i 

^ Jo Jo Jo 



aia 2 . . . a n _i(l - {a% + a 2 + . . . + a„-i)) 



(aia 2 . . . a„_i + (1 — (01 + a 2 + • • • + a n -±))(a 2 a3 . . . a„_i + a 3 a 4 . . . a„_iai + . . . + aia 2 . . . a„_ 2 )) 



2 



1 



2 n-3 



dxi / dx 2 . . . / dx„_ 2 
Jo Jo 

2 



1 + (xi + X 2 + ■ ■ ■ + Xn-z) 2 1 + X\ + X 2 + ■ ■ ■ + X n - 2 

: — log 1 

2{X\ + X 2 + ■■ ■+ X n - 2 ) 1-{X 1 +X 2 + ..- + X n - 2 ) 

iXi... X n -2 (1 - (Xi + X 2 + ■ ■ ■ + In-2) 2 ) 

4(xi + . . . + X n - 2 )xi . ..X n -2 + (1 - (xi + ... + X n ^ 2 ) 2 )(x 2 X 3 . . .X n -2 + ••• + X 1 X 2 ■ ■ . X„_ 3 ) 



(65) 



Finally if in (164p and (|65p one changes notations Xi — > a% one has shown 

/>00 />0O />00 />0O />00 

2™~ 1 / du uKq(u) u = J dai • ■ • / da n _ 2 / da„_i / da n — e~ 

1 



70 JO JO JO 

1 rl — Ui pl—Un-3 pl—U n -2 

dai / da 2 . . . / da„_ 2 / da„_i 

Jo Jo Jo 



W n -i + (1 - u n -i)v n -\ 
(66) 

1 + U n -2 4 



1 1*1—141 /*1 — tt n — 3 

dai / da 2 . . . / da„_ 2 log 

Jo Jo 1 - u n-2 4u„_ 2 W„_ 2 + (1 - U„_ 2 jU„-2 

(67) 



and 



w 



2" 3 / d««%«= / dai... / da„_ 2 / da„_! / da„ — e 

JO JO JO JO JO V n 

1 pl—Ul rl—U n — 't />1 — Un—2 /1 \ 

, f A f A f A W n -i{l-U n -l) 

dai / da 2 . . . / da„_ 2 / da„_i — — 

K-l + (1 - Un-i)V n -x) 



Jo Jo Jo 



(68) 



dai / da 2 . . . / da„_ 2 — log 1 



Jo Jo 



2u„_ 2 1 - u„_ 2 / (4u„_ 2 u>„_ 2 + (1 - u 2 n _ 2 )v n - 2 y 



(jnU) and indicate that J °° du uKo{u) n and J °° du u 3 ifo('") TI are periods. This 
is the case in general for / °° du u m Ko(u) n with m odd. 



5 CONCLUSION : WHY LOOKING AT BESSEL 
INTEGRALS ? IS £(5) IRRATIONAL ? 

In Section @ some arguments (and examples) were given of why weight k Bessel 
integrals are expected to be related to weight k — 1 (and below) zeta numbers. 
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One has found the surprisingly simple PSLQ identity 

C(5) =PSLQ J o U K 0( U ) S dU -^J " 3 K °^ 8 dU ^ 

If one could prove for k — 8 the irrationality claim on (|28p . in particular on 
J °° u Kq(u) 8 du and v? Ko(u) s du, and derive ((69|) . then one would have a proof 
of the irrationality of £(5). A related issue is to generalize the algebraic construction 
of Section to weights k > 5, with Bessel basis of dimension > 3 (see Appendix 
B). 

Coming back in Section (|3.3[) to C(3), C(2) and ipi(l/3) — -0i(2/3), it would be 
rewarding to find a systematics beyond the cases listed and have simple expressions 
(Apery- like numbers) for the y ra 's satisfying the recurrence relations associated to 
these continuous fractions with initial conditions {y(0) = 1, 2/(1) = 0} or {y(0) = 
0, y(l) = 1}. With the Apery initial conditions {y(0) = 1, y(0) = q}, the recurrence 
y(k + 1) - (2k + l)(3fc 2 + 3fc + l)y(k) + k 6 y(k - 1) = associated to the continuous 
fraction 

z{k - 1)= (2k + im^+3k + i) + Z (k) z(0) =PSLQ m 1 

has for solution [TU] 

fc! 3 * /^A 2 / 2l s 2 



y(A){„ ( o)=i,»(i)=i} = ^fc2^ ) [ k J ( 7 °) 

where fc! 3 factorises. Similarly for the k = 4 continuous fraction 

— 16fc 6 12 
Z(fc " 1} = (2fc + l)(5fc 2 + 5fc + 2) + z(fc) Z(0) = W) " 2 (71) 

fc! 3 A /fc\ 2 /2i\ (2{k-i) 



(see in [5] the solutions for the second ((2) and the k = 3 continuous fractions listed 
in Section (13.31) ) . 

One would also like to integrate u K$(u) n du further up to possibly one integra- 
tion left on an integrand which should contain a term of the type (log(l + x)/(l — 
x)) n ~ 2 times a rational function of x yet to be determined, possibly allowing for a 
derivation of Q69p. 

Finally the fact that Bessel integrals fall in the category of periods might also be 
an indication of a deeper meaning yet to be understood. 

Acknowledgments: S.O. acknowledges some useful conversations with S. Mashke- 
vich in particular for helping in the numerics involved in the PSLQ searches of 
Appendix C. We also would like to thank Alain Comtet for a careful reading of the 
manuscript. 
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APPENDIX A: INTEGRATING FURTHER 



1) Trying to integrate J °° du uI (u)K (u) n 
Using again the integral representation (JS3J) and (it = 2y/i) 



one can integrate over t using 



to obtain 



4-k 

*>(«)=£ w (73) 

k=0 y '' 



dt t k e~t* = * (74) 



Introducing as above the variable j3, integrating over a n and then trivially over 
one finally obtains 

/>oo 

T- 1 / duuI (u)K (u) n 
Jo 

dai / da 2 ■ ■ ■ da n - 2 / da n 



1 W„_ilt„_i + (1 - U n -ijV n - 

(75) 



V 1 1 / du «/oW^o(«)" = / da x ... I da„_ 2 / da„_i / da„ — e u « 



a result to be compared to ([571) . 



One can push the integration one step further following the same procedure as in 
Section to obtain an expression again in terms of u n - 2 ,v n ^ 2 and tu n -2 but 
somehow more involved than (|67[) since its denominator contains a square root. 



2) More on integrating J °° du uK (u) n , J °° dw M 3 i^o(w) n and J °° du uI (u)K (u) 
n = 3 



dwwi4T (u) 3 = 



1 

22 



dsi log 



1 + asi 



l-si 4(a?i)xi + (1 - (a;i) 2 ) 



= / dxilog 



1 + x 1 1_ 

1 — x\ 1 + 3x 2 



du uKq(u) 



1 

2^ 



dxi 



l — Xl 



dx 2 log 



1 + X± + X 2 1 



1 - X\ - X 2 X\ + X 2 AX\X 2 + (1 - (Xl + X 2 ) 2 ) 

(76) 
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n = 5 



I duuK (u) 5 =— / dxi / dx 2 / dx 3 log 
/o ^ Jo Jo Jo 



1 - (Xi + £ 2 + X 3 ) 



&{x\ + x 2 + x 3 )x\X 2 x 3 + (1 - (xi + x 2 + x 3 ) 2 )(x 2 x 3 + x 3 xi + Xix 2 ) 

(77) 



Integrating further 
n = 3 



3 

o Jo 



1 + X\ 1 f 1 1 + Xi 



duuK {u) = / dxilog- = / dacilog- 

l-Xil + 3^ 7 l-xi (l+a;i) 3 (l + (^-) 3 ) 



1 f 1 , 1 + u . 

2 A du TT^ logu 



-/ d«(l +u)^(-l) fc M 3fc l0g'W 
2 J t_ n 

Vi(l/3)-Vi(2/3) 

(78) 



fc=0 

2^ V ; V (3/c + l) 2 (3/c + 2) 2/ 12 



where one has made the change of variable u = \ + x x x x and used 

1 ni 



d««*(Iogfi)" = (-l)*^^ 



n = 4 

1 f 1 , Z" 1 ^ 1 , , I + Z1+Z2 1 



y dwuA'o(w) 4 = ^ y dxi y dx 2 log 



1 — X\ — X2 Xl +X2 1 — {x\ + X 2 ) 2 + 4X1X2 



i/2 



= — - f dx— log — - — / dw 

2 3 2 7 x & 1- X J_ X y l-y 

= - / da;- log / dy 1 

2 jo 1 1 - ^ Jo 1 - y 

11 Z 1 , l n l + x, 2 7C(3) 
where one has made the change of variables x = x\ + x 2 and y = x\ — x 2 and used 

^rC s =2^— + -i — = o / dx:c ^^+^vr =7 — Try / —(log, — s 1 

2 s 1 rt s n s 2 / n — 1 + 1 (s — 11! 7 n x 1 — x 

n—l 



23 



or rather, from 

f°° &uuK (uf = - f dx-(log^) 2 = i f du-^(. 
Jo 4 J x 1 - x 2 J l-u z 



\ / 1 d u (log U ) 2 £ u 2 « = ^ (81) 

•*° n=0 



with the change of variable u — 



l+x • 

oo 



One knows for example that duu 3 Ko(u) A is a linear combination with rational 
coefficients of 1 and duu 3 Ko(u) 4 namely 

/>oo />oo 

4/ duuK Q (u) 4 - 16 / duii^oM 4 = 3 (82) 
Jo Jo 



One has 



A 3^/,4 /\ /" 1_X1 , ^ 1 + (^1+X 2 ) 2 ] 1 + Xl +X 2 \ 1-{X 1+ X 2 f 

&UU 6 K (u) = dX! dx 2 —, — - — r— log- r - 1 

Jo Jo V 2{x 1 +x 2 j l-(xi+x 2 ) J 



(xi + x 2 ) 



j (83) 



1 — (xi + x 2 ) 2 + Ax\X 2 



so that 



\ r 1 -* 1 , . i+ xl + X2 i 2 

dxi / dx 2 log ■ 



JQ 



1 — x\ — x 2 x\ + X2 1 — (xi + x 2 ) 2 + 4xiX 2 



. 1 + (xi + x 2 ) 2 1 +xi +x 2 _ \ 1 - (xi + x 2 ) 2 4xi x 2 „ 

2(X!+X 2 ) ° g l-(x!+X 2 ) y (X!+X 2 ) 2 / X 2 ^ 



1 - (xi + x 2 ) 2 + 4xix 2 
(84) 

has to be satisfied. With the change of variables x = x\ + x 2 and y = x\ — x 2 

f 1 , Z -1 -* 1 , A + (xi+x 2 ) 2 , 1 + Xi+X 2 A 1 - (xi +x 2 ) 2 4xix 2 
7o 7o ^ { 2(x 1+ x 2 ) l0g l-(x 1+ x 2 ) " 1 J (x 1+ x 2 ) 2 



1 - (xi + x 2 ) 2 + 4xix 2 



2 f 1 dx (±±* log i±£ - 1^ r d„ +y){x - y) 



2J ~ V 2x "° 1 - x 'J x 2 y w (1 - y 2 ) 2 

dx {-&- bg — - l ) — 2 + (1 + a } 2 bg —x 

\ 1-x 2 ^ 1 + x 2 l + x \ 2 

^^x-^^x- 108 !-^- 1 ; (85) 

so that (|54"|) becomes 

^ j A/i 1 + ^2 ,1-X 2 ,1 + X 2 , 1 + X , 2 \ „ , o „. 

/ dx -log- ) 2 -4 — log- 1) J =3 (86) 

/ n V x 1 — x x 2x 1 — x / 
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which is indeed true. 
Finally 

i y /"l /"I— cti 2 

duuIo(u)Ko(u) 4 — — dai / da2 



n 



4 7o Jo oia 2 (ai +a 2 ) — (oi +C2)(1 — ai - 02) 

(87) 

With the change of variables x = a% + 02 and y = a\ — 0,2 one gets 



&uuIq{u)Kq(u) =— / da; / dy 



1 Z" 1 , , 4 1 1 



4 Jo Jo xx— 2 + yx — 2 — y 
Ax— r / dy- 



2 7 x(x-2) J a x-2 + y x-2-y 
If 1 . 1 



1 r 1 _L 3C(2) 

- / du -log it = 

2.L 1-u 2 8 8 



where u = 1 — x. 



3) K-Bessel integrals summary 

/•OO 

/ duui^otw) = 1 
JO 







1 AT 

duJCo(u) = | = -^(Vo(l/3)-Vo(2/3)) 



1 

duuK a (u) 2 = - 







(i 



3 

d W X„( U ) 2 = -C(2) 



Jo 12 

duK {u) 3 =? 

/•OO y 

y duuK Q (uf = -C(3) 
diiA" (w) 4 =? 
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? 



du uK {uf -72 du u'K^uf = PSL Q 77C(5) 
o Jo 



on — 1 roc 

lim / duuK (u) n = e 2 ^ 



n— foe Ti 







lim — / du/r (u) n = 2e*« 

n->oo n! Jq 

where V'o(l) is minus the Euler constant. 



APPENDIX B: HIGHER ORDER RECURSIONS 

From the iteration (flTJ)) it is easy to get for k > 5 higher order recursions generalizing 
the recursions (|4"2l l?T|) for k = 4, 3. For example in the k — 4 case, (|16[) or (|3"3")l 
becomes, defining x(k) — I^ Q , y(k) = 1^. 4 , 

x(fc) ^ _ l + 2fc f d(fc) ^ f x(fc + 1) ^ 



y(fc) ; (fc- i)fc(fe + 1)2 V a(fc) c(fc)/^(Hl) 

where a(k) , b(k) , c(k) , d(k) are given in (|3"5)) . 

One can for example start from the inverse iteration 

/ x(k + l) \ _ ( a(k) 0(h) \ ( x(k) 
\ y(k + l) J { 7 (fc) S(k) J { y(k) 

where 



(90) 



(fc-l)fc(fc + l) 2 1 
W 1 + 2* b{k)c{k) - a{k)d(k) C[ ' 

(fc-l)fc(fc + l) 2 1 dW 

PW ~ l + 2fc b(k)c(k) - a{k)d(k) { ' 

(k~l)k(k + l) 2 1 
71 j f + 2fc b{k)c(k) - a(k)d(k) a[ ' 

6[k) - l + 2k b(k)c(k) - a{k)d(k) m (Ji) 

Consider k ->• fc + 1 in ([501) to get 

x(fc + 2) = a(k + l)x(k + 1) + £(* + l)y(k + 1) (92) 
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Invert (JHOJ) 

*(*) = Q(fc) , (fc) ! /3(fc)7(fc) (^)^ fc + 1) - mv(k + D) (93) 

and use ([93[) to eliminate y in ([9"2"| . One obtains a recursion for x 

/3(fc)a;(fc+2)-(/3(fc)a(fc+l)+/3(fc+l)J(fc))x(fc+l) + (a(fc)(5(fc)-^(fc)7(fc))/3(fc+l)x(fc) =0 

(94) 

that is to say, using (f9"Tj) 

fc 4 x(fc - 1) - (2fc - l)(2fc + 1)(2 + 5fc + 5k 2 )x(k) + 16(1 + fc) 2 x(fc + 1) = (95) 

When multiplied by 16fc 2 the recursion ([55]) is (g2|) for .T(fc) = 8 fe (2fc)!fc!x(fc). 

In the n = 5 case, the iteration (|T^|) for x(fc) = I^koi y(k) = ^2fc 2 an< ^ z (^) = -^2^4 
is 

5(l+2fc) 10(l+2fc) n 

2(fe+l) (fc+1) U 

l+2fc (l+2fc)(9+17fc) 3(l + fc)(l+2fc) 

k Ik 2 k 2 

n 6(l + fc)(l+2fc) (l + fc)(l+2fc)(-8+13fc) 

U fc(fc-l) 2fe(fe-l) 2 

Again one may start from the inverse iteration 

a(k) p{k) 7 (fc) 











y(fcH 






V z(fcH 


-i) j 





7j(fc) 0(fc) l(fc) 




(97) 



where in particular f3(k)((k) — j(k)e(k) = and 8(k)9(k) — e(k)r)(k) = to account 
for the vanishing elements of (|M|) . In ([57)) consider fc — > fc + 1 

z(fc + 2) = r](k + l)x(k + 1) + 9(k + l)y(k + 1) + t(fc + l)z(k + 1) (98) 

y(k + 2) = <y(fc + l)x(fc + 1) + e(ft + l)y(k + 1) + C(A + + 1) (99) 
and k -> fc + 2 

2/(fc + 3) = <5(fc + 2)a;(/c + 2) + e(fc + 2)y(/c + 2) + ((k + 2)z(k + 2) (100) 
([551 [TUP]) allow to eliminate the z's in ([9"B"]l 

+ 3) - 5(fc + 2)x(fe + 2) - e(fc + 2)y(fc + 2)) 



C(k + 2Y 

= n(k + l)x(k + 1) + 6(k + l)y(k + 1) + ffi + (y(fc + 2) - *(* + l)x(fc + 1) - e(fc + l)y(fc + 1)) 

C(k + 1) 

(101) 



Invert (197 



= (-c(fc)g(fc) + e(k)L(k)) x(k + 1) + ( 7 (fc)g(fc) - mm± y(k + 1) f1(121 

1 ' -f(k)5(k)8(k)-a(k)((k)9(k)-p(k)6(k)L(k)+a(k)e(k)i(k) [ ' 
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and use (|102p to eliminate the y's in (|10ip . One finally obtains a recursion for x, 
which reads, using the matrix elements in f|96[) 



8k 5 x{k - 1) - (-1 + 2k) ^4 (5 + 28fc + 63k 2 + 70k 3 + 35k 4 ) x(k) 

- (1 + fc)(l + 2k) (2 (285 + 518fc + 259fc 2 ) x{k + 1) - 225(2 + fc)(3 + 2k)x(k + 2)) ^ = 

(103) 

Generalizing this construction to higher order recursions for k > 5 is straightforward 
-see [TT] for recursions of this type (see also [12] for an example involving £(5) and 
C(3)). 

APPENDIX C: A FAMILY OF DOUBLE NESTED 
INTEGRALS 

Coming back to (|2I9[) use 

poo poo pu 

g(u)du / f(x)dx = / f(u)du / g(x)dx 
a Ju Jo Jo 

to rewrite l p i a s as 

/>oo />U 

I p 2 a e=8 duuKaiufiuK^u)] 2 dx xK 1 {x)I 1 {x)K {x) 2 
Jo Jo 

-4 / duuK^ufK^u) 2 I dxxK (x)xK 1 (x)(xK 1 (x)I (x) - xK {x)h(x)) 
Jo Jo 



+ / uKo(uf{uKx(u)Ydu (104) 
Jo 

Next use xii (a:)if (x) + xIq(x)Ki(x) = 1 so that 

/•OO f-U 

/ P 2 a6 =8/ duu^K^ufKxiuf \ dxxK 1 (x)I 1 (x)K (x) 2 



Jo JO 

/>oo />u 

+8 / duuKoiufK^u) 2 / dxx 3 ^!^)/!^)^^)^ 
Jo Jo 

-4 / duuA' (u) 2 /fi(w) 2 / dxa; 2 ^^)^!^) 
Jo Jo 

poo 



+ / uKoiuriuK^ydu (105) 



o 
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One has for the next to last term in (1105 



du«^(ti) 2 ]fi()i) 2 / dxx 2 K (x)K 1 (x) = 2 / u K {u) 2 K^u) 2 ((uK^u)) 2 - l)du 

2 r°° 



wA'i(u) ((uKxiu)) 2 - l)dA (u) 3 
K (u) 3 d(uA 1 ( W ))(( U X 1 ( W )) 2 - 1) 



uK (u) 4 du-2 / uK (u) 4 (wAi(m)) ^dw 
(106) 



so that 



7 p 2 q6 =8/ duw d KoiuyK^uY I dx iX 1 (i)J 1 (a;)if (a;) i 
Jo Jo 

duuA' (w) 2 A'i(ii) 2 / da;a; 3 Ai(j:)/i(.T)Ao(.T) 2 
Jo 

ziAo(u) 4 dw 

uA (w) 4 (uAi(u)) 2 dii (107) 



2 

3 7 



Defining 

f n (x) = x n K {x) 2 K 1 {x) 2 and ffn (a;) = x n K (x) 2 K 1 (x)h(x) 
equations flU [T07|) imply that 

C(/ 3 ,.9i) + C(fi,93) =pslq^ J uK (u) 6 du-—J^ u 3 K (u) 6 du 

_ JL £( 3 ) _ J!L £( 5 ) (108) 
96 sv ; 1280 sv ; v ; 

as confirmed by a direct PSLQ check. 

So the meaning of © when rewritten as (|108[) might be that it is the symmetric 
form CCfsiSi) + C(/iiff3) which is a linear combination with rational coefficients 
of simple Bessel integrals of weigth 6 and £ numbers like £(5) and below. This 
suggests to look at other symmetric sums of double nested integrals sharing this 
property. A PSLQ search confirms that C(/5; <?i) + C(/ij5s) belongs indeed to this 
category 

211 r „, , 6 ,„. , 3953 r,3^ 6 . 



C(/5,3i) + C(/i,55) =pslq YY^Tj J ^o(«) b d u + y u d ifo(«) b d« 

111 93 

+ C(3)-— — C(5) (109) 

9216 9 ^ v ; 5120 sv ' K ' 

as well as C(/ 7 ,fli) + C(/i>07) 

«*•»'> + W'-^=^« msooo /„ » A '»("' 6d « + iS5eooo/ « 3A '»<«) 5d » 

+ lls-|« 3 >-I> < 110 » 
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and Cl/3,55) + C(/5,53) 



C(/3,5 5 ) + C(/5,53) =pslq — ?8921_ f°° rf( ) 6 d 1151533 f 00 3 6 
su^MiTsus.w 691200 i vy 1382400 J v; 

14653 25 A/n , 279 

+ 184320 + 192 20480 C(5) < m > 

We infer that for any two odd positive integers n, to one should have that £(/„, ffm) + 
C(/m, <?n) can be rewritten as a linear combination of 1, J Q uKo(u) 6 du, J Q u 3 Kq(u) 6 
and £(3) and C(5). Finding the coefficients of the linear combination can be in prin- 
ciple obtained by generalizing integration by part procedures in [TJ[21[11[5] to double 
nested integrals. 
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